Intrinsic Electric Dipole Moments of Paramagnetic Atoms: Rubidium and Cesium 
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The electric dipole moment (EDM) of paramagnetic atoms is sensitive to the intrinsic EDM con- 
tribution from that of its constituent electrons and a scalar-pseudo-scalar (S-PS) electron-nucleus 
interactions. The electron EDM and the S-PS EDM contribution to atomic EDM scales as ~ Z A . 
Thus, the heavy paramagnetic atomic systems will exhibit large enhancement factors. However, the 
nature of the coupling is so small that it becomes an interest of high precision atomic experiments. 
In this work we have computed the EDM enhancement factors of the ground states of Rb and Cs 
due to both the electron EDM and the S-PS EDM using the relativistic coupled-cluster (RCC) 
theory. The importance of obtaining the precise enhancement factors and the experimental results 
in deducing a reliable limit on the electron EDM is emphasized. 

PACS numbers: 11.30.Er, 31.15.bw, 31.30.jp 



In the contemporary era of Large Hadron Collider 
(LHC) in which physicists worldwide are waiting with 
excitement and hoping to find answers to many out- 
standing questions haunting the field of particle physics 
for decades, non-accelerator probes like the observation 
of the electric dipole moment (EDM) of any elementary 
particle attempting to look for new physics is quite sig- 
nificant. It is indeed complementary to the efforts of 
LHC, however, at the low energy regime. In fact, despite 
the relentless search for a non-zero EDM of elementary 
particle like the electron for more than 50 years, no con- 
clusive results have been obtained so far. This is quite 
intriguing and a challenging task for the experimental- 
ists because of the nature of the couplings involved. The 
EDM of the electron if observed would be extremely small 
(< 1.6 x 10~ 27 e — cm [lj) and hence it certainly falls into 
the category of the highest precision atomic physics ex- 
periments. The extensive search for EDMs has resulted 
in the application of several ingenious experimental tech- 
niques including a variety of cooling and trapping tech- 
niques, electric and magnetic shielding mechanisms etc. 
These latest low energy high precision EDM experiments 
on Rb and Cs atoms are being pursued in different labo- 
ratories 0, H, |i[ and many proposals are still being con- 
sidered including an experiment to be conducted in the 
microgravity environments on-board space missions 
All these proposed EDM experiments are striving hard 
to improve the sensitivities by 2-3 orders of magnitude 
higher than the current limits. It would be possible to 
obtain better limit on the electron EDM (e-EDM) by 
combining the precise theoretical and experimental re- 
sults. In this context, we have carried out high precision 
numerical calculations of the EDM enhancement factors, 
first in the series, of Rb and Cs which are sensitive to 
the intrinsic EDM of the electron and the scalar-pseudo- 
scalar (S-PS) electron- nucleus interactions 

The e-EDM has far reaching consequences ranging 



from cosmology [7] to particle physics Q. The non-zero 
observation of EDM in any non-degenerate system will be 
a direct proof of time-reversal violation in nature. Var- 
ious models of particle physics have predicted different 
values for the magnitude of EDM of the electron, many 
of which are close to the current experimental limits un- 
like the value predicted by the Standard Model which 
cannot be achieved even in a foreseeable future and thus 
the observation of non-zero EDM can be helpful in con- 
straining different models of particle physics which pre- 
dict the size and nature of CP violation which, in turn, 
has wider implications such as understanding the baryon 
asymmetry of the universe masses of certain heavy 
particles, of course, depending on the type of the model 
considered Thus EDMs undoubtedly will unfold a 
novel direction for understanding the physics beyond the 
Standard Model. 

We employ the open-shell relativistic coupled-cluster 
(RCC) theory with single, double and a subset of leading 
triple excitations in our calculation which is known by 
the name CCSD(T) method in the literature. The gen- 
eral discussion of the CCSD(T) method is described in 
[To| and the references therein. The method employed in 
particular in calculating the EDM enhancement factors 
is discussed in detail in JjJ, ll2| . However, for the sake of 
completeness, we briefly outline the procedure here. 

We first obtain a reference wavefunction (|$o)) for the 
closed-shell state with N — 1 electrons, N is the total 
number of electrons in the system including a valence 
electron (v), by solving the Dirac-Fock (DF) equations; 
the Hamiltonian of which is given by, 

+ (fa - l)mc 2 + Vv(rO} + ]T Vfc(r«)(l) 

i i<j 

where c is the speed of light in vacuum, a and (3 are the 
Dirac matrices, Vj^f is the nuclear potential and Vc is the 
Coulomb interaction term. 



2 



In the framework of the RCC theory, we construct the 
exact wavefunction (|>Ei> )) for the corresponding valence 
electron system as, 



(2) 



where is the excitation operator for the core (occu- 
pied) orbital electrons, corresponds to the excitation 
operator for the valence and valence-core electrons and 
|$„) = ajjl^o) is the single valence reference state where 
a\ is the particle creation operator. 

Even in the absence of any external field, any param- 
agnetic atom may have two dominant sources of EDMs; 
one, the intrinsic e-EDM contribution and the other due 
to P- and T- violating S-PS electron-nucleus interactions. 
They are neglected due to their meager influences on de- 
termining any physical quantities. However, their effects 
are considerably larger in heavy atomic systems and can 
be studied to probe some subtle effects. The intrinsic 
e-EDM contribution is purely quantum field theoretical 
in nature and will be non-vanishing only in the relativis- 
tic treatment of the problem [l8} . By considering these 
interactions into account, the total atomic Hamiltonian 
can be written as, 



H = H 



H 



EDM 



(3) 



where Hq is the unperturbed Hamiltonian given by 
Eq.([T]) and Hedm is the part of the Hamiltonian per- 
turbed either by intrinsic e-EDM contribution or by the 
S-PS EDM contribution and they are given respectively 
by, 



H EDM = - " X! 



and 



(4) 



(5) 



where d e is the coupling constant for e-EDM, 7 s is the 
product of the Dirac matrices i.e, 7 s = ij 1 ^ 2 ^ 3 ^ , pj 



is the momentum vector of the j th electron, Gf is the 
Fermi coupling constant, Cs is the dimensionless S-PS 
constant, A is the mass number of the atom and p_\f is the 
nuclear density. The Hedm mixes the atomic states of 
opposite parities but with the same angular momentum. 
As its strength is sufficiently weak, we consider only up 
to the first-order perturbation in wavefunction. Thus, 
the modified atomic wavefunction for the valence electron 
state V is given by, 



(6) 



where A = d e for e-EDM, whereas, A = GfCsA) \[2 for 
S-PS EDM. 

In the RCC ansatz, the cluster operators for calculat- 
ing the perturbed wavefunctions are given by, 



T = T(°) + A T«; and S v = + A 



(7) 



where and Sv are the first-order corrections to 

the unperturbed cluster operators and , respec- 
tively. In the relativistic CCSD approximation, the per- 
turbed and unperturbed cluster operators are taken to 
be, T = Ti + T 2 and S v = S V1 + S V2 . We have excited 
all the core electrons to all possible virtual states in the 
present calculations. To obtain better accuracies of the 
wavefunctions we have calculated the perturbed and un- 
perturbed cluster amplitudes with a convergence of less 
than 10 -6 in all the cases. 

The expectation value of an atomic EDM (D a ) for the 
state \^f v ) is given by, 



D„ 



(*v\D\Vv) 

(*i 0) |^ 0) ) 



(8) 



where D = ef is the electric dipole moment operator. 
The final expression for the EDM enhancement factor 



Do 
dc 



or S 



D„ 



= ) in terms of the perturbed 



G F C s A/y/2> 

and unperturbed cluster amplitudes is given by, 



R{ovS) 



{{1 + S^}DW{TV> + T«^ 0) + + {Sl iy * + S^T^ + T« 1 }£("){! + S< 0) }} |*„) 



<$„|e 



T(»)' e T(°) 



where £K°) 



P (0)t 



"(0) 



D e .In the above expression the 



DF term is contained in £)C°) Si. It corresponds to con- 
sidering D in and one order of EDM interaction and 
zero order of residual Coulomb interaction in Sv . It can 



I 

be explicitly written in the DF approximation as, 
(*v\Hedm\$i)(*i\S\$ v ) 



(9) 



(10) 
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where the intermediate states differs from the va- 
lence reference state |<&„) which is 5s for Rb and 6s for 
Cs either by a single valence or core orbital and of oppo- 
site parity. 

It has to be noted that the accuracy of the coupled- 
cluster results depend on the size of the configuration 
space of the RCC wavef unctions considered. However, an 
increase in the configuration space leads to a larger mem- 
ory requirement and also computational time. Thus, one 
is limited by the maximum size of the configuration space 
that one has considered in an actual calculation. Our 
preliminary results with a medium sized basis set and 
with a limited number of multipoles satisfying Coulomb 
interaction selection rules were published before in 11|. 
However, here we have taken more than 5 different sets 
of large configuration spaces with all possible multipoles 
and systematically studied the nature of the correlation 
trends of the dominant RCC terms in Eq.fTU]) such as, 

LWS^, D@)TW, and S| 0)t P)^ 1 ' with their hermi- 
tian conjugates contributing to the e-EDM enhancement 
factor (R) for the Rb atom and plotted the variation in 
the magnitude of the correlation contribution of these 
terms as (a), (b) and (c) in Figure 1, respectively. The 
RCC formulation contains all-order relativistic Coulomb 
correlation effects like core-correlation through the term 
£)(o) TW, pair-correlation and core-polarization through 



FIG. 1: The variation in the trends of the dominant RCC 
terms such as, (a) DW (b) DWT W (c) Sf )+ flMsi 1 ' 
and (d) Total e-EDM enhancement factor (R) with the size 
of the basis set for Rb. 
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the term £K°) . The difference between DS£ } and TABLE I: The EDM enhancement factors of the ground state 



the DF terms is the largest of all the correlation contri- 
butions. It is clear from all the above three diagrams that 
the magnitude of correlation increases with the increase 
in the size of the basis set and more or less gets satu- 
rated beyond 95. It is interesting to observe the same 
convergence trend in the e-EDM enhancement factor R 
shown in the last diagram (d) with the size of the basis 
set. Thus, we obtain the convergence in the magnitude 
of EDM enhancement factor by taking sufficiently large 
number of basis functions to represent the space spanned 
by the RCC wavefunctions. Using the same basis sets we 
also obtain the EDM enhancement factors due to S-PS 
electron-nucleus interactions (S) for Rb. We have com- 
puted the EDM enhancement factor R for a few basis 
sets of size larger than 100 for Cs. The value of S for 
Cs is taken from [l^] • Finally the converged results both 
for Rb and Cs along with the Dirac-Coulomb results are 
tabulated in Table I. We have also compared our results 
with some of the published ab-initio methods based on 
finite order MBPT [HQ, MCDF+MBPT [HG3 and 
semi-empirical one-electron theories 17 1. 

The EDM matrix element is sensitive to the wavefunc- 
tion near the nuclear region where as the dipole matrix 
element is sensitive to the wavefunction at far-off regions 
from the nucleus. Thus there is a greater need for more 
reliable relativistic calculations of the atomic wavefunc- 
tions of both the nearby and far-off states, from the nu- 
cleus. The RCC theory employed in the present work 



of Rb and Cs atoms. R (= 



, , and S (= „ T^n-js ) are the 
enhancement factors due to electron EDM and scalar-pseudo- 
scalar electron-nucleus EDM respectively. 



Atom 



R 



Reference 



Rb 



Cs 



DF 


Total 


DF 


Total 




19.55 


25.74 


11.91 


15.64 


This work 




24.6 
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16.47 
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94.19 


120.53 






This work 
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72.44 
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114.9 






[13] 




114 






[14] 




130.5 






[16] 




119 
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has such a potential. The accuracy of EDM results will 
depend on; (i) the excitation energies, as the reliable val- 
ues of these are crucial as the EDM enhancement factor 
inversely depends on the energy differences as shown in 
Ea. tfTOl) . (ii) the spontaneous electric dipole (El) transi- 
tion amplitude for the transition from the lowest 2 P\/2 
state to the ground state and (m) the magnetic dipole hy- 
perfine constants (A) of the lowest 2 S , 1 / 2 and 2 P\/2 states. 
To a good approximation, the error in the EDM matrix 
elements can be taken as f 



A2 S ■ As P . We have 
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TABLE II: The excitation energy (in cm 1 ) and the El tran- 
sition matrix element (in a.u. ) for 5 2 S —* 5 2 Pi/2 transition 
and the magnetic dipole hyperfine constants of the 5 2 Si/2 
and 5 2 Pi/ 2 states of Rb atom ( 85 Rb; r = 5/2"). 



Property 


Transition/ 


Mag 


mitude 




State 


This work 


Expt. [Ref.J 


Excitation energy 


5 2 S _ +5 2p i/2 


12579.87 


12578.95 [19] 


El tr. amplitude 


5 2 S^5 2 P 1/2 


4.26 


4.23 [20] 


Hyperfine Constant 


5 2 S 1/2 


1009.33 


1011.91 [21] 


Hyperfine Constant 


5 2 Pi/ 2 


117.71 


120.64 [22] 



hancement factors of their corresponding ground states 
and comparing them with the accurately known experi- 
mental results. 

We acknowledge Dr. Raj at Kumar Chaudhuri of Indian 
Institute of Astrophysics, Bangalore, India for develop- 
ing certain parts of the codes. We also thank the staff 
of CD AC Terra-flopp Super-computing Facility at Banga- 
lore, India for their support and the computational facil- 
ity for carrying out this work. 



shown all the above discussed property results for Rb in 
Table II and compared them with those of the published 
experimental results some of which are known up to fairly 
high accuracies such as, for example, Rb ground state 
hyperfine energies are considered for the frequency stan- 
dards and our results agree reasonably well with them. 

The differences between the computed and the accu- 
rately known experimental results presented in Table II, 
give the individual errors and by adding the errors in 
quadrature we estimate that the maximum error in our 
EDM results could be around 0.7%. We have also ob- 
tained the Cs results of similar quality. 

The limit on the c-EDM will be obtained by combin- 
ing the precise values of both the atomic EDM measure- 
ments and the theoretical enhancement factors. The bet- 
ter e-EDM limit from Cs EDM experiments has resulted 
from [23| . By combining their experimental result of Cs 
EDM and our computed enhancement factor we obtain 
the limit (—1.5 ± 3.2) x 10~ 26 e — cm. However, the best 
ever e-EDM limit is obtained from Tl EDM experiment 
as (6. 9 ±7.4) x l(T 28 e- cm by Regan et al. Q. There are 
a few state-of-art EDM experiments which are currently 
underway to measure atomic EDMs in Rb and Cs 
aiming to achieve almost 2-3 orders of higher sensitivity 
than the reach of the current experimental setups and 
when these attempts come to fruition, by using our en- 
hancement factors a better limit on intrinsic e-EDM can 
be achieved. 

In conclusion, we have calculated the EDM enhance- 
ment factors of the ground states of the paramagnetic 
atoms such as, Rb and Cs atoms due to both the intrin- 
sic electron EDM and the scalar-pseudo-scalar electron- 
nucleus interactions using the RCC theory which inher- 
ently has an all-order relativistic many-body nature. We 
have systematically studied different correlation terms, 
which contribute dominantly to atomic EDM, for a num- 
ber of basis sets of varying configuration spaces. We have 
estimated the conservative upper limit of the error bar 
in our EDM results to be 1% by studying various other 
atomic properties of the lowest 2 «S , 1 / 2 and 2 P\/2 states of 
Rb and Cs which dominantly contribute to the EDM en- 
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